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Abstract 

Despite ablation and drag processes associated with atmospheric entry of 
meteoroids were a subject of intensive study over the last century, little atten¬ 
tion was devoted to interpret the observed fireball terminal height. This is a 
key parameter because it not only depends on the initial mass, but also on the 
bulk physical properties of the meteoroids and hence on their ability to ablate 
in the atmosphere. In this work we have developed a new approach that is 
tested using the fireball terminal heights observed by the Meteorite Observa¬ 
tion and Recovery Project operated in Canada between 1970 — 1985 (hereafter 
referred as MORP). We then compare them to the calculation made. Our 
results clearly show that the new methodology is able to forecast the degree of 
deepening of meteoroids in the Earth’s atmosphere. Then, this approach has 
important applications in predicting the impact hazard from cm- to meter¬ 
sized bodies that are represented, in part, in the MORP bolide list. 
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1 Introduction 


Deriving the meteoroid trajectories in the atmosphere is of particular interest 
to researchers. On the one hand, orbital parameters can be derived based on the 
time of appearance, meteor position and initial velocity, allowing us to estimate their 
parental relationship with parent bodies - asteroids and comets. On the other hand, 
a knowledge of physical parameters such as mass, velocity, deceleration, height, etc, 
at different points of its trajectory turns out to be very useful so as to predict the 
energy of a possible surface impact, locate meteorite fall and/or understand the 
ablation and other mass loss mechanisms occurring along the flight. 

Various photographic and video techniques have been developed to obtain the 
most accurate and systematic observations of meteors. Whipple and Jacchie (1957) 


(later modihed by 

McCrosky and Posen 

(1968) 

Pecina and Ceplecha 

(1983 

) and 

Ce- 

plecha et ah 

(1993 

)) derived a methodology by considering the separate meteor trails 


obtained when shuttering the video image. This technique allowed them to study the 
problem at shorter flight intervals. This analysis was dependent on the body prop¬ 
erties’ average values provided by the bibliography. However, a reliable theoretical 
flight model is still to be developed. The accuracy of theoretical results usually re¬ 


quires a very good precision in observation techniques. For example, Ceplecha et ah 


(1993) developed a theoretical model which included meteor fragmentation, but it 
needs very precise hreball records. 

Along the years different methods have been developed in order to increase the 
accuracy of the theoretical models. 

One of the hrst theoretical models for meteors (known as Single Body Theory) 
was developed by Hoppe (1937). It was an extensive study of the flight mechanics and 
thermodynamics processes. Levin (1956, 1961) studied the meteoroids’ atmospheric 
entry with account for fragmentation and deceleration. He concluded that the mass 
of the body is related to its middle section by means of a parameter that characterizes 
the rotation of the hreball. 

Later on, Ceplecha and McCrosky (1976) explored which hreballs of the Prairie 
Network were ordinary chondrites. The authors considered the hreball terminal 
height as a main characteristic factor. It was concluded that carbonaceous material 
shall ablate more readily and, consequently, these bodies may have shorter trajec¬ 
tories. The authors derived an empirical criterion (equation ([^) that established a 
weighted relation between the hreball terminal height and other hight properties of 
the hreball (namely air density at terminal height, pe [s'/cm^], the preatmospheric 
mass, moo [d], preatmospheric velocity, V/o [km/s], and the zenith distance of the 
meteor radiant, Zr [degrees]). 
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PE = log pE + A log moo + B log Voo + C log {cos Z^) 


( 1 ) 


This expression also has a theoretical meaning based on the single body theory, as 
explained in their work. Coefficients A, B and C are obtained by using a least squares 


£t to 156 fireballs of the Prairie Network (McCrosky and Boeschenstein 1965). Owing 


to this criterion (equation (p(|)), [Ceplecha and McCrosky (1976) classified the Prairie 
Network fireballs into four different groups. In their discussion they suggested that 
ordinary chondrites should all belong to the same range of PE values (—4.60 < PE). 


Besides, Ceplecha and McCrosky (1976) took advantage of the fireballs’ ob¬ 


served properties in order to shed more light in the validity of the previous cri¬ 
terion (equation ([^). Two parameters were used: K (the shape-density coefficient) 
[cm^^f^^/^jand a (ablation coefficient) The average values (for all the ob¬ 

servational measurements) are used to define a new parameter, SD: 


SD = < log K > + < log a > 


( 2 ) 


The average numbers, < log K > and < log a >, are weighted by the ratio of 
the deceleration to its formal rms error. Then, SD is a parameter that has little 
influence from observational errors. As Ceplecha and McCrosky (1976) stated, SD 


depends on the second derivative of the observational measurements via these two 
refined numbers, whereas PE is chosen as the simplest possible empirical expression. 
Therefore, SD is a new criterion which can be compared against the PE criterion. 
These two parameters ([^ and ([^ turned out to be related when the meteor initial 
mass could be considered small or the ablation was large. 


Slightly different methodology was suggested by Wetherill and Revelle 


They included four meteorite selection criteria to build up their classification. 


(1981). 


Wether 


ill and Revelle (1981), gave more importance to the dynamic mass than to the pho¬ 


tometric mass; besides, they also took into account the deceleration of the body and 
the light curves to identify the survived meteorites among the fireballs registered by 
the Prairie Network. The authors highlighted the importance of the observed ter¬ 
minal height concluding that for meteorite-producing fireballs its value should agree 
with the theoretical value, calculated using dynamic mass, as well as with that of 
Lost City fireball to within 1.5 km, when scaled for mass, velocity, and entry angle 
in accordance with classical single body meteor theory. 


In Revelle (1979), the study of the interaction between large meteoroids and the 


atmosphere is done via a quasi-simple ablation model. The results are compared 
to photographically recorded meteorite falls as well. Later publications have gath¬ 
ered and expanded the physical problem of the deceleration of meteoroids in the 


atmosphere (e.g. Bronshten (1983)). 
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Halliday et al. ( 1989a||b ) studied observed fireball properties to derive the pres¬ 
ence of correlations. They used 44 MORP recorded hreballs to classify as strong, 
moderate, weak or not having any correlation the observed data (i.e. initial velocity, 
total light emitted by the hreball, initial and end masses, initial and end heights, 
orbital elements, etc.). Despite of the observational errors (cameras not able to him 
all the trajectory, not clear sky, etc.), they found some strong correlations: mass 
lost by ablation versus the peak brightness, and duration of luminosity recorded by 
MORP versus zenith distance of the radiant. 


Stulov et al. (1995), Stulov (1997) and Gritsevich (2007) proposed a new method¬ 


ology. Instead of using the average values as input parameters, they gather all the 
unknown values into two new variables a (ballistic coefficient) and (3 (mass loss 
parameter), mathematically introducing similar idea with scaling of parameters as 
suggested by |Wetherill and Revelle (1981). Adjusting the resulting equation to the 
trajectory observed, these new variables can be derived for each meteoroid. The 
resulting values allow to describe in details the meteoroid trajectory in the atmo¬ 


sphere and invent new classihcation scale for possible impacts (Gritsevich et al. 2011 


Gritsevich et al. 2012). This allows to determine other important parameters, such 


as preatmospheric and terminal mass values, ablation and shape change coefficients, 
as well as terminal height. The methodology to determine terminal height has been 
implemented for fully ablated fireballs by Gritsevich and Popelenskaya (2008). In 


the present study we significantly specify and expand the applicability range of this 
methodology by testing it on larger data set. 

In the following sections we present the results of applying this last methodology 
to a large number of MORP hreballs, including suspected meteorite-producing events 
included in the table 6 by Halliday et al.| ( 1996| ). Alternatively, we suggest a more 
accurate method of calculation. Section takes a look on previous and present 
terminal height determinations. We then compare observed values to our derived 
values in section Section contains our discussion. Finally the conclusions and 
suggestions for future research are presented in section 


2 Theory 

The equations of motion for a meteoroid entering the atmosphere projected onto 
the tangent and to the normal to the trajectory are well known. Once we consider 
some simplihcations 


see 


Gritsevich (2010)) we have: 


= --CdPaV^S 

dt 2 ^ 


( 3 ) 
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—Vsin'f 


( 4 ) 


dh 
dt 

M is the body mass, V is its velocity, t is the time, h is the height above the 
planetary surface, 7 is the local angle between the trajectory and the horizon, S is 
the area of the middle section of the body, pa is the atmospheric density and q is 
the drag coefficient. 

Equations ([^, Q are complemented by the equation for the variable mass of the 
body: 


H 


.dM 

dt 


= -^ChpaV^S, 


( 5 ) 


where H* is the effective destruction enthalpy, and Ch is the heat exchange coef¬ 
ficient. Two other expressions are added to solve the problem. On the one hand, we 
suggest that the atmosphere is isothermal, and so p = exp{—h/ho), where ho is the 


scale height. On the other hand, from the research carried out by Levin (1956| 1961) 


we can assume that the middle section and the mass of the body are connected by 
the following relation s = m^, where p is a constant. The parameter p characterizes 
the possible role of rotation during the flight and can be calculated based on the ob¬ 


served brightness of a hreball using the method proposed by Gritsevich and Koschny 


(2011). In Gritsevich (2008c) different values for p are discussed. If p = 0 there is 


no body rotation, whereas if p = 2/3, the ablation of the body due to its rotation is 
uniform over the surface, and the shape factor does not change. Generally we have 


0 < p < 2/3. According to the recent results (Bouquet et ah 2014) the majority of 
the MORP fireballs were found to have p values closer to 2/3. 

We now introduce dimensionless variables, M = Mpm, V = VeV, h = hoy. 
Pa = Pop, S = SeS. Here, the subscript e indicates the parameters at the entry to 
the atmosphere, ho is a planetary scale height (we use 7.16 ■ 10 ^ m for the Earth) 
and po is the atmospheric density near the planetary surface. 

Taking into account all these previous considerations, we can turn the equation 
system ([^-(|^-([^ into the following equations (see Gritsevich 2009): 


dv 

m— = apvs 
dy 

dm ^ 9 

—— = 2afdpv s 
dy 


( 6 ) 

( 7 ) 


These equations have been expressed as a function of two variables: 
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( 8 ) 


And 



PohoSe 

MeSin'j 


/? = (1 -/i) 


2cdH* 


(9) 


a is called the ballistic coefficient. It is proportional to the mass of the atmo¬ 
spheric column with the cross section along the trajectory divided by the body 
preatmospheric mass; it characterizes the drag intensity. 

(3 is known as the mass loss parameter. It is proportional to the fraction of kinetic 
energy of the unit mass of the body that is transferred to the body in the form of 
heat divided by the effective destruction enthalpy. 

The analytical solution of the equation system (|^-Q with the initial conditions 
y = oo, V = 1 and m = 1 is: 


m = exp[—{l — — /i)] 


( 10 ) 


y = ln2a + (3 — InA, A = Ei{f3) — Ei{(3v‘^), 


( 11 ) 


Ei{x) = 


e^dt 


( 12 ) 


For quick meteors, one can neglect deceleration and set the velocity as constant 
for the luminous segment of the trajectory, then v = V/Ve = 1. A strong evaporation 
process takes place so (3 becomes high (/3 >> 1). This means that the meteoroid 


ablation occurs rapidly in the absence of drag. For this situation, Stulov (1998, 2004) 
developed an asymptotic solution: 


n = 1, ^ = 1 — 2a(3e ln2a(3 < y < oo 


(13) 


It should be noted that relation (13) does not describe the drag process until the 
point where m = 0. The meteor velocity begins to decrease in a certain vicinity of 
this point. In order to account for this change in velocity it was proposed to combine 


equation (13) suitable for high 13 values with equation (10) suitable for arbitrary f3 


values (see Gritsevich (2008c) and references therein). Such substitution of equation 


(10) into equation (13) fully escapes dependency of y and gives the following height 


dependence of velocity at high j3 values: 
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V = 


ln{l — 2a(3e 


+ 1 


1/2 


In2al3 < y < oo 


(14) 


It is notable, that since condition of changing velocity was applied in the deriva¬ 
tion of eqnation (14), the equation ( [l4) ) is no longer applicable for constant-velocity 
meteors. Thus if condition v = V/Ve = 1 is strictly set then solely equation (13) 
should be considered. 

In the opposite case, when /3 = 0, the analytical solution (@-(0 becomes 
(|Gritsevich||2008cP : 


m = 1, y = Ina — ln{—lnv) 


(15) 


2.1 Brief explanation of a and /3 derivation 


In their work, Gritsevich (2007) explained the best way to derive a and (3 pa¬ 
rameters. Herein, we summarize the main points of the methodology. For a deeper 
explanation of the whole methodology we recommend a careful reviewing of the 
mentioned articles. 


As explained in Gritsevich (2008b) a weighted least-squared method is applied 


to the observed height and velocity values using equation Assigning manually 
the weighted factors may be quite complicated, so, since the height and velocity of a 
meteor decrease while it gets closer to the Earth’s surface, the solution was proved 
to perform better if we take an exponential form of equation ( [IT| . Mathematically 
it means that instead of direct comparison of yi (observed) with yi (calculated) the 
comparison is made between values. The obtained results are then compared 
with the original data and a very good agreement is shown (see e.g. 

(2008c)). Then, equation (pT| turns into the following form: 


Gritsevich 


2aexp{—y) — Aexp{—/3) = 0, 


A = Ei{^) - Ei{/3v^), 
-, , r Edt 


( 16 ) 

(17) 


a and 13 are derived by minimizing the following expression (18). We take into 
consideration that the analytical trajectory expression obtained shall accurately ad¬ 
just the list of observed trajectory {yi,Vi) values. 
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Q{a,P) = '^{Fi{yi,vi,a,p)f 


( 18 ) 


2=1 


where 


Fi{yi, Via, 13) = 2aexp{-yi) - Aiexp{-I3) (19) 

In equation (18) we assume n > 3, i.e. there are at least 3 points along the 


trajectory available for the calculation, one of which should contain initial velocity 
value, 14 , needed for scaling. It is also assumed that deceleration was noticeable 
during the luminous trail, i.e. there are at least 2 other than beginning points 
{k,j < n) along the trajectory such as Vk < Vj < I 4 . In principle, recent proliferation 
of instruments and methods dedicated to observing meteors makes it possible to fulhll 


this condition even for such small meteor particles as the Draconids (Borovicka et al 


2007 


i Kero et al. 

2012 

Ye et al. 

2013; 


Vaubaillon et al. 2015). In ordinary cases. 


when equally reliable data are provided along the observed trajectory (e.g. Halliday 


et al. (1996)), we suggest the use of all available {hi, Vi) values, but this is not a strict 
requirement and specihc initial hlters on observational data may apply depending on 
studied case. 


3 Results 

We have derived terminal heights by means of the previous development. In 
order to do so, we have used the a and [3 values which have been previously derived 


using the methodology described above by Gritsevich (2009) for the MORP hreballs. 


These heights are compared to the hreballs observed terminal height values; |Halliday 
et al.| (1996) collected them along with other hreball parameters such as: brightness, 
beginning and ending velocity, beginning height above sea, etc. 

The results of the following calculations are shown in table The hrst row of 
this table contains the ID name of the hreball. Column 2 contains I 4 which is the 
velocity of the body when it enters the atmosphere. Column 3 shows the terminal 
velocity observed, Vj. Values for Vj and Vj are taken from Halliday et al. (1996). 


The a and [3 values appear in columns 4 and 5 respectively, and they are taken from 
Gritsevich (2009). The values of hj, hn and hjjj in columns 6, 7 and 8 correspond 


to the terminal heights derived in this study as we will see in the following lines. 


Finally, in column 9 the observed terminal heights from Halliday et al. (1996) are 
listed for comparison with the values derived in our study. 




































Table 1: Terminal Heights 


Fireball No. 

14 , km/s 

14 , km/s 

a 

/3 

hj, km 

hi I, km 

hill, km 

hobs, km 

18 

18.5 

5.7 

24.13 

1.475 

30.54 

32.73 

29.65 

27.6 

123 W 

16.3 

8.7 

37.22 

1.111 

31.61 

35.92 

33.29 

32.6 

138 

16.9 

9.5 

38.90 

2.889 

38.77 

39.84 

37.83 

37.0 

141 

32.1 

29.0 

707.19 

6.403 

65.24 

67.87 

67.28 

65.1 

144 

23.5 

13.5 

137.13 

1.301 

42.08 

45.96 

43.54 

42.9 

169 

22.9 

9.3 

50.52 

1.575 

36.30 

38.54 

35.71 

34.0 

172 [“1 

12.5 

8.1 

13.13 

3.379 

32.12 

33.20 

31.52 

31.2 

177 

15.8 

8.7 

77.63 

1.777 

40.24 

42.69 

40.33 

39.8 

187 

18.4 

14.0 

83.76 

3.815 

46.25 

47.86 

46.57 

46.3 

189 

14.5 

4.8 

34.47 

0.757 

28.32 

33.42 

30.02 

27.9 

192 

21.0 

14.6 

94.01 

2.458 

43.93 

46.29 

44.59 

44.4 

195 [“1 

25.2 

7.6 

35.22 

1.486 

33.30 

35.45 

32.37 

30.4 

204 H 

13.0 

8.7 

10.13 

3.620 

30.75 

31.80 

30.21 

29.5 

205 [“1 

19.7 

7.5 

36.84 

0.716 

28.39 

33.99 

30.70 

28.9 

207 W 

17.9 

6.5 

24.45 

0.775 

26.03 

31.14 

27.83 

25.9 

218 

18.5 

13.2 

64.65 

0.339 

27.07 

40.49 

38.53 

37.9 

219 [“1 

18.4 

7.8 

12.51 

2.060 

28.23 

29.69 

27.11 

26.1 

223 [“1 

27.1 

9.5 

18.33 

1.809 

30.03 

31.67 

28.84 

27.1 

225 

21.2 

10.0 

75.92 

0.706 

33.47 

39.64 

36.64 

35.4 

229 

12.3 

6.1 

43.73 

4.564 

42.88 

43.12 

41.46 

40.7 

231 [“1 

27.9 

11.7 

52.72 

1.371 

35.61 

38.41 

35.52 

34.2 

232 

35.0 

29.8 

434.10 

0 



56.58 

56.5 

235 

19.1 

17.5 

171.01 

14.343 

60.85 

61.60 

61.19 

60.7 

241 

16.7 

9.8 

50.27 

0.878 

32.08 

38.00 

35.51 

34.9 

245 H 

13.5 

8.9 

83.60 

0.478 

31.37 

41.68 

39.44 

39.0 

261 [“1 

12.4 

7.8 

39.13 

1.084 

31.80 

37.04 

34.79 

34.4 

268 [“1 

19.1 

15.7 

39.44 

5.157 

43.02 

44.51 

43.54 

42.9 

276 [“1 

23.5 

6.5 

18.15 

1.237 

27.24 

29.99 

26.74 

24.4 

285 [“’^1 

14.54 

2.7 

8.25 

1.700 

23.87 

25.41 

22.30 

19.8 

287 

23.4 

10.5 

40.85 

1.750 

35.53 

37.57 

34.91 

34.6 

288 W 

12.4 

4.1 

9.05 

1.215 

22.13 

25.09 

21.93 

20.2 

299 [“1 

23.6 

14.7 

48.36 

1.042 

33.03 

38.41 

36.12 

35.7 

300 

21.5 

10.2 

202.36 

2.071 

48.20 

49.80 

47.33 

46.2 

301 

36.1 

24.4 

693.65 

0.662 

48.85 

57.42 

55.30 

55.6 

303 [“1 

14.1 

8.2 

44.15 

1.260 

33.74 

37.82 

35.41 

34.9 

304 

16.2 

6.6 

48.17 

2.442 

39.10 

40.10 

37.64 

36.8 


Continued on next page 


9 



Table 1: Terminal Heights (continued) 


Fireball No. 

14 , km/s 

14 , km/s 

a 

/3 

hj, km 

hi I, km 

hill, km 

hobs, km 

307 ^“1 

21.0 

3.8 

12.08 

1.760 

26.85 

28.29 

25.21 

22.0 

310 

17.0 

10.5 

67.74 

0.559 

30.98 

39.79 

37.36 

36.8 

313 [“1 

16.7 

10.2 

134.20 

0.432 

34.03 

44.33 

41.83 

41.1 

314 

14.6 

9.3 

159.93 

0.118 

26.00 

45.27 

42.83 

42.1 

331 

13.3 

7.0 

37.94 

0.598 

27.32 

34.81 

31.98 

31.1 

340 

29.1 

25.6 

1165.35 

5.894 

68.22 

70.41 

69.70 

69.1 

345 H 

17.4 

10.6 

46.21 

0.829 

31.07 

37.52 

35.11 

34.3 

346 

15.7 

8.5 

21.81 

1.997 

31.98 

33.99 

31.67 

31.0 

364 H 

11.3 

5.3 

20.34 

0.567 

22.47 

29.84 

26.77 

25.4 

384 

21.1 

12.1 

476.52 

2.121 

54.50 

56.47 

54.28 

53.8 

388 

21.9 

14.8 

81.67 

0.542 

32.10 

41.88 

39.74 

39.4 

390 

19.7 

16.7 

958.47 

1.344 

56.24 

64.51 

63.44 

63.4 

391 

19.5 

9.0 

72.46 

1.323 

37.63 

40.75 

37.96 

36.9 

410 

31.0 

27.3 

901.17 

2.912 

61.33 

66.59 

65.78 

64.9 

414 

29.9 

27.7 

375.26 

12.707 

65.61 

66.90 

66.49 

66.5 

419 

24.7 

16.2 

329.97 

0.962 

46.21 

52.38 

50.23 

49.8 

425 

17.6 

12.3 

43.41 

1.699 

35.76 

39.65 

37.83 

37.2 

436 

29.0 

25.9 

224.06 

3.797 

53.27 

57.73 

57.02 

57.0 

445 

21.3 

13.7 

119.22 

0.700 

36.64 

44.44 

42.16 

41.6 

462 

19.1 

10.6 

140.03 

2.015 

45.36 

47.40 

45.13 

44.5 

481 

13.7 

5.9 

77.70 

1.046 

36.45 

40.43 

37.44 

35.8 

486 

33.2 

26.4 

624.60 

1.059 

51.46 

59.57 

58.16 

58.1 

503 W 

14.8 

10.5 

25.27 

1.765 

32.15 

36.01 

34.25 

34.1 

511 W 

18.1 

9.8 

48.38 

1.156 

33.78 

37.95 

35.36 

31.7 

516 

20.5 

13.3 

104.87 

3.703 

47.65 

48.54 

46.92 

46.3 

518 

15.5 

8.7 

75.32 

2.391 

42.15 

43.70 

41.55 

40.8 

521 

18.8 

12.8 

188.55 

0.998 

42.46 

48.77 

46.74 

46.4 

528 

19.5 

8.4 

37.25 

1.525 

33.89 

36.33 

33.54 

32.4 

529 

18.5 

10.2 

35.26 

1.779 

34.60 

37.05 

34.69 

34.0 

536 

12.9 

9.3 

61.78 

3.332 

43.11 

44.72 

43.24 

42.8 

549 

26.5 

12.1 

111.65 

0.740 

36.57 

42.39 

39.35 

38.1 

555 

25.5 

24.1 

481.92 

11.668 

66.79 

69.22 

68.88 

68.8 

557 

36.5 

27.8 

138.83 

2.139 

45.73 

49.47 

48.00 

47.9 

559 

26.0 

20.1 

135.61 

3.454 

48.99 

51.04 

49.76 

49.6 

561 

16.0 

11.7 

154.14 

2.817 

48.45 

50.70 

49.18 

48.6 

565 [“1 

13.2 

10.7 

12.57 

4.913 

34.49 

35.95 

34.93 

34.8 

567 [“1 

23.4 

8.1 

101.02 

1.560 

41.19 

43.28 

40.32 

38.7 


Continued on next page 
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Table 1: Terminal Heights (continued) 


Fireball No. 

14 , km/s 

14 , km/s 

a 

/3 

hj, km 

hi I, km 

hill, km 

hobs, km 

577 

22.0 

21.6 

86.30 

33.029 

61.92 

64.52 

64.40 

64.3 

589 

24.3 

12.0 

167.95 

1.549 

44.78 

47.44 

44.83 

43.9 

593 

20.3 

16.1 

313.35 

1.224 

47.56 

54.78 

53.38 

53.1 

598 

17.9 

11.1 

118.26 

1.471 

41.90 

45.61 

43.41 

42.9 

610 

29.7 

25.7 

384.29 

1.053 

47.94 

58.39 

57.42 

57.4 

615 

20.5 

8.8 

55.62 

1.434 

36.32 

38.98 

36.15 

34.9 

626 [“1 

13.50 

10.0 

31.59 

0.118 

14.38 

35.57 

33.74 

33.5 

628 

28.9 

21.6 

600.34 

1.850 

55.17 

59.35 

57.78 

57.6 

635 

35.5 

18.5 

96.63 

1.883 

42.22 

44.32 

41.90 

41.2 

654 

13.8 

6.9 

43.15 

1.356 

34.10 

37.31 

34.65 

33.8 

661 

30.9 

17.5 

218.48 

1.281 

45.30 

49.19 

46.74 

46.1 

664 

22.1 

13.0 

145.75 

2.257 

46.46 

48.32 

46.21 

45.8 

669 W 

20.6 

8.6 

28.80 

1.434 

31.60 

34.22 

31.36 

30.6 

672 [“1 

13.7 

10.0 

22.79 

1.243 

28.90 

34.77 

33.03 

32.7 

683 [“1 

17.6 

9.3 

37.46 

1.340 

33.00 

36.43 

33.86 

33.1 

687 [“1 

16.7 

5.9 

42.83 

0.534 

27.37 

34.43 

30.97 

28.9 

702 

28.7 

19.9 

306.92 

0 



48.52 

48.2 

708 

15.7 

10.8 

63.51 

1.727 

38.60 

42.29 

40.42 

40.1 

715 

31.7 

27.9 

141.98 

6.700 

54.06 

55.85 

55.16 

54.8 

725 

29.5 

28.0 

93.72 

14.748 

56.74 

58.63 

58.33 

57.9 

727 

31.3 

27.7 

1160.73 

2.100 

60.80 

68.00 

67.20 

67.2 

741 

31.9 

28.5 

370.61 

3.764 

56.81 

61.32 

60.61 

60.6 

751 

18.9 

8.9 

44.19 

0.377 

25.10 

34.91 

31.78 

30.5 

752 

13.6 

11.9 

8.67 

9.418 

36.48 

37.32 

36.69 

36.5 

763 

27.6 

25.2 

254.30 

4.410 

55.24 

59.93 

59.34 

59.4 

769 

17.4 

10.5 

50.00 

1.233 

34.47 

38.84 

36.52 

36.0 

771 [“1 

17.1 

7.8 

41.75 

0.731 

29.44 

35.32 

32.28 

31.0 

774 

33.0 

27.1 

85.66 

0.994 

36.78 

45.99 

44.74 

44.7 

776 

13.1 

10.0 

33.21 

5.497 

42.25 

43.01 

41.89 

41.6 

792 H 

17.3 

7.7 

37.43 

0.388 

24.12 

33.57 

30.34 

28.9 

819 

17.7 

11.1 

123.90 

1.034 

39.71 

45.18 

42.90 

42.4 

823 

25.0 

19.4 

32.61 

4.927 

41.33 

42.42 

41.28 

40.0 

829 

28.9 

15.1 

359.02 

0.869 

46.08 

51.51 

48.76 

47.7 

835 

24.2 

8.2 

66.13 

2.474 

41.46 

42.31 

39.77 

38.4 

840 

23.6 

6.1 

22.35 

1.766 

31.28 

32.81 

29.81 

27.7 

843 

36.5 

33.7 

642.92 

4.526 

62.07 

67.22 

66.69 

66.6 

844 H 

14.9 

9.8 

23.89 

1.694 

31.46 

34.91 

32.92 

32.6 


Continued on next page 
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Table 1: Terminal Heights (continued) 


Fireball No. 

Ve, km/s 

Vt, km/s 

a 

/3 

hj, km 

hji, km 

hni, km 

hobs, km 

852 

15.9 

4.6 

21.27 

2.541 

33.53 

34.26 

31.69 

29.9 

857 

18.5 

12.0 

62.52 

1.224 

36.02 

40.87 

38.74 

38.3 

860 

23.3 

17.6 

330.46 

0.283 

37.46 

52.98 

51.25 

51.1 

861 

13.3 

7.2 

49.82 

2.720 

40.11 

41.24 

39.14 

38.6 

865 

26.5 

18.1 

258.85 

0.564 

40.65 

50.29 

48.20 

47.8 

868 

27.8 

25.6 

1234.31 

9.416 

71.99 

73.94 

73.47 

73.4 

871 

23.9 

20.2 

134.32 

0 



47.87 

47.8 

872 H 

14.8 

5.8 

12.34 

0.631 

19.66 

25.98 

22.68 

20.8 

873 

16.5 

14.8 

927.96 

5.170 

65.65 

68.89 

68.24 

68.0 

878 

15.5 

14.1 

58.70 

1.512 

37.08 

47.62 

46.96 

46.9 

883 

26.1 

8.8 

95.53 

2.115 

42.97 

44.16 

41.45 

40.1 

884 H 

12.4 

8.3 

26.37 

2.569 

35.15 

37.13 

35.35 

34.9 

886 

13.8 

7.4 

55.66 

0 



33.18 

32.0 

888 [“1 

25.5 

9.0 

31.85 

1.171 

30.87 

34.06 

30.92 

29.2 

892 

28.9 

6.2 

32.77 

1.514 

32.92 

34.84 

31.65 

28.9 

901 

32.6 

27.3 

636.32 

0.422 

45.01 

60.28 

59.10 

59.0 

902 

29.6 

22.4 

447.06 

2.547 

55.35 

58.29 

56.85 

56.7 

920 

16.7 

12.7 

318.85 

4.052 

56.26 

57.69 

56.42 

56.3 

921 

28.1 

22.4 

1077.86 

8.721 

70.47 

70.77 

69.96 

69.3 

925 [“1 

26.4 

4.2 

42.87 

1.503 

34.79 

36.67 

33.42 

29.8 

929 

35.0 

30.1 

342.27 

2.944 

54.48 

58.95 

58.03 

58.0 

930 H 

13.4 

9.5 

38.75 

1.052 

31.51 

37.94 

36.06 

35.6 

933 

15.4 

11.1 

43.80 

1.779 

36.15 

40.12 

38.41 

38.0 

937 [“1 

17.9 

6.4 

40.19 

0.881 

30.50 

34.96 

31.69 

29.8 

947 

24.9 

24.0 

771.47 

32.917 

77.58 

78.31 

78.13 

77.9 

966 

31.3 

27.8 

1384.16 

0.604 

53.14 

68.34 

67.51 

67.5 

967 

30.0 

15.0 

239.38 

0 



43.09 

41.9 

977 [“1 

23.3 

9.4 

30.90 

1.065 

29.98 

33.76 

30.69 

29.4 

984 

36.2 

26.1 

128.24 

0 



43.00 

42.7 

995 

28.9 

26.0 

240.98 

4.700 

55.31 

59.07 

58.42 

58.8 

996 

26.9 

14.1 

59.56 

2.133 

39.65 

41.36 

39.04 

38.4 

1010 

18.4 

10.1 

71.58 

1.378 

37.84 

41.28 

38.79 

38.1 

1011 

37.6 

36.6 

546.56 

24.155 

72.90 

75.26 

75.10 

75.0 


[“1 Probable meteorite fall according to Halliday et al. (1996). This calculation corre¬ 
sponds to the Innisfree meteorite and it is based on the data from Halliday et al. (1981), 


see 


Gritsevich (2008c). 
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We have followed a step by step methodology. In hrst instance, terminal height 
is derived as the height corresponding to the trajectory point where m = 0 based on 


equation (13). The substitution of condition m = 0 into equation (13) leads to: 


hi = ho ■ yt = ho ■ ln{2a(3) 


( 20 ) 


where ho = 7.16 ■ 10^ m is the scale height, and yt is the scaled terminal height. 
As it can be seen in table differences between observational and analytical 
values are remarkable (the standard deviation is a = 4.11 km) and we conclude that 


the method proposed by Gritsevich and Popelenskaya (2008) cannot be extended 


to be applicable for decelerated hreballs without further modihcation, though it is 


usable to derive parameters of small meteoroids with constant velocity (Popelenskaya 


2010[ ). Figure]^ shows graphically the obtained differences between the observed and 
derived terminal height values based on equation (20). As a reference, we have also 


plotted the line for which hj = hobs {hots are the terminal heights observed taken 


from Halliday et al. (1996)). We then took into consideration that v shall vary up 


along the trajectory. Terminal heights were obtained from the equation (14): 


hii = ho-yt = ho- In 


2al3 


( 21 ) 


(l - 

Once again, results can be found in column 7 in table Notice that Vt = Vt/Vf. 
(where Vt is the terminal velocity observed, and W is the velocity of the body when 
it enters the atmosphere) is the scaled terminal velocity. As we can see, though 
differences between derived and observed terminal heights (h//, hobs) have decreased 
(the results’ standard deviation is a = 1.52 km) there is still a wide spread in 
the results. However, if we represent terminal heights derived vs terminal heights 
observed (£gure[^, a global linear tendency can be found. This indicates that this 
methodology seems to obtain a reasonable first-order approximation. 

One possible explanation for the quantitative differences obtained in these results 


is the use of the approximated function (14) established for high jS values instead 


of using more general equation For easy comparison with equation @ , the 

equation (14) can be rewritten as: 


y = \n 


2al3 


( 22 ) 


(l_e/3h^-i)) 

We have investigated either any improvement to this approximation itself, i.e. 


any possibility of adjusting the equation (14) to better match equation (11) could 


also bring better alternative solution to the problem studied here. To answer this 
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question, we have used results recently obtained in (Gritsevich et al. 2015). The 
authors have examined the possibility of improving approximated function (|14[). The 
functions ( pT| and ( [2^ do not have singularities and they are monotonous on the 
interval 0 < u < 1. The dependency on a is the same (additive as ln{a)) for 
both, equations ([IT|) and (22). Under these conditions a shift along parameter (3 


could quantitatively improve approximation represented by equation (22). And, in 


particularly, based on further functional analysis, it was found that much better 
approximation for the function ([IT| can be achieved with the use of the function: 


y = \n 


2a (/9- 1.1) 


(23) 


(1 _ e(/9-i.i)(^2_i)^ 

On hgure ([^ we have graphically demonstrated the quality of the approximation 
relied on equation (23). We compare functions |/(/d, v) = y{a, (3, v)—ln{2a), where y is 
taken according to equations (0, (0 and ( [2^ . The comparison made for different 
(3 values conhrms that approximation (23) has better accuracy than previously used 


equation (22). In order to improve the results of the present study, we derived a new 


expression to calculate the terminal height of a hreball based on the referred above 


approximation (Eq. (23): 


hill = ho ■ yt = ho ■ In 


2a(/3- 1.1) 


(24) 


_ g(/3-l.l)(«2-l)j 

A hrst justihcation to take up this further step is supported by hgure it rep¬ 
resents the differences {hi — hobs) vs. (3. As expected, equation (20) seems to give 
better results for larger f3 values. The same happens with equation (21), though in 
this case this dependency to higher (3 values is smaller (as can be noticed from its 
standard deviation a = 1.52 km). 

Nonetheless, this new modihcation turns out to provide higher accuracy to our 
results. Now, the standard deviation becomes a = 0.75 km. We can notice it 
graphically in hgure As we can see, this hgure shows the total agreement between 
derived and observed terminal heights using this new formula. Besides, it even allows 
us to derive terminal heights in those few cases where f3 values are null. 


4 Discussion 


As explained in section]^ equation (20) should lead to good results when (3 
values are high. In hgure we notice that for (3 > 5 the diherences between derived 
and observed terminal heights are small. However, this equation did not take into 
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hi 


Figure 1: Representation of derived terminal heights using equation (20), hi 
the observed terminal height, hobs- Line indicating hi = hobs is also plotted. 


versus 
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hii 


Figure 2: Representation of derived terminal heights using equation (21), hn versus 
the observed terminal height, hobs- Line indicating hn = hobs is also plotted. 
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/3 = 2 


/3 = 4 



t; = \ yv; 


/3 = 6 



CN 

I 


3 

2.5 
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1.5 
1 

0.5 
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r = vyi' 

3= 10 

.^ 


l_I_I_I_I 
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1.5 - 
1 - 
0.5 - 

0 -^^^^- 1 

0 0.2 0.4 0.6 0.8 1 

i' = vy v; 


Equation (11) 
Equation (22) 
Equation (23) 


Figure 3: Graphical comparison of functions 0. (0 and for the selected 
3 values. The authors note, that for real hreballs, terminal velocity Vt = Vt/Ve is 
usually higher than the one shown on the graphs (see experimental Vt and 14 values 
given in table [^, so only range < u < 1 would be applicable for practical purposes. 
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Figure 4: Representation of the differences between derived and observed terminal 
heights {hi — hobs) versus [5 values. 
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hill 


Figure 5: Representation of derived terminal heights using equation (24), hm versus 
the observed terminal height, hobs- Line indicating hm = hobs is also plotted. 
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account the decrease in velocity close to the terminal point of the trajectory. Given 
that the amount of hreballs studied is large, we continued our analysis considering 
this fact (see equation 21). 

Formula (21) seemed to provide accurate results with regard to the observations; 
in fact, a clear global lineal tendency was found on the graph with the observed 
and calculated values of terminal height. Nonetheless, remarkable discrepancies are 
present along the whole range of (3 values. No exceptions for higher, lower or inter¬ 
mediate /3 values (see table [^. 

The modihcation made by introducing new recently proposed approximation 
(Gritsevich et al. 2015) in our methodology leads to higher accuracy. This result 
may create some astonishment. Given that observations could carry different kind 
of errors: bad astrometric measurements, incomplete recording of the luminous tra¬ 
jectory of the meteor, blurring effects, etc., it should be expected that some of the 
analytical values differed from the observed ones. This could be explained if we 
remind that the differences obtained are of the order of 1 km. Most of the afore¬ 
mentioned observational errors shall introduce smaller distance errors so these effects 
shall be comprised in our calculations. 

Additional observation can be made based on figure ([^ that equation (24) pro¬ 
duces its best £t for hjjj for the fireballs which terminated at higher altitudes. This 
result is expected since fireballs with survived terminal mass penetrate lower in the 
atmosphere, and they have, on average, lower values of mass loss parameter (3. For 
example, Halliday et al. (1989a) derived dynamic mass estimation only late in the 


trail for those MORP hreballs that penetrated below a height of 40 km. Gritsevich 


(2008a) calculated dynamic terminal masses of MORP hreballs based on the whole 
observed trajectories regardless of the terminal height value; indeed, all larger ob¬ 
jects and/or suspected meteorite droppers penetrated below 40 km, with an average 
terminal height of 29 km. Thus, for planetary defense applications one can note, 
that larger objects could initiate hreballs penetrating deeper in the atmosphere than 
our methodology would predict. 

It is also notable, that meteor height may be expressed as a function of time. 
Thus, the ability of predicting a terminal height may be directly linked with the 
forecast of a total duration of meteor phase. This can create scope to a new class 
of problems, such as, for example, insights into determination of luminous efficiency 
based on meteor duration and calculation of critical kinetic energy needed to produce 
luminosity. In future, several other tightly related problems may be considered, 
including a possibility of initial meteor height determination and deriving parameters 
a and (3 without the beginning part of a meteor trail. Furthermore, the existing 
solution to the inverse problem of deriving these parameters described in section 
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of this study is only applicable for the decelerated hreballs. There are still large 
amount of archived data and many new meteor registrations for which the recorded 
deceleration was not signihcant and/or it is in the order of observational errors. For 
these cases observed terminal heights will be crucial in setting margins on possible 
a and /? values based on equation (20). Some initial steps in solving this problem 
were taken in Popelenskaya| (|2010[). 


Moreover, even for obviously decelerated and well-studied cases, such as Innisfree 
meteorite (hreball No. 285 in table [^, published terminal heights may differ depend¬ 
ing on a taken data-reduction approach (e.g. 21 km in Halliday et ah (1981); 19.8 


km in Halliday et ah ( 1996[ )) and, therefore, any additional theoretical constrains on 
the terminal height and hreball duration may be very helpful. 

Finally, it is interesting to compare our results with PE criterion introduced by 


Ceplecha and McCrosky (1976) reviewed in the Introduction section of our study. In 


our methodology we used a different approach which, in particular, avoids A, B, C 
coefficients needed to be empirically set in PE criterion. We operated with dimen¬ 
sionless parameters a and [3 which by their dehnition (equations ([^ and (|^) gather 
all the unknowns. The terminal height was a key sought parameter without passing 
to air density calculated at terminal height. We adjusted a theoretical expression 
which can be used to calculate terminal height. To verify the proposed dependency 
of terminal height on parameters a and f3 we calculated terminal heights for MORP 


hreballs with known (observed and reported by Halliday et al. (1996)) terminal height 
and known values of a and {3 parameters. However, if one would substitute a and 
f3 in our criteria with their dehning expressions (equations ([^ and ([^), one could 
hnd many similarities with PE criterion. Indeed, a depends on the preatmospheric 
cross-section-to-mass ratio (the ratio which can be easy converted into dependency 
on bulk density, preatmospheric mass and shape coefficient) and trajectory slope 7 
related to Zj^. Mass loss parameter is proportional to preatmospheric velocity in 
power of 2 and inversely proportional to the effective destruction enthalpy. Thus, 
surpassing unknowns with described in our study methodology we hnd similar set of 
physical parameters ahecting the degree of deepening of meteoroids in the Earth’s 
atmosphere. 


5 Conclusions 

Previous work has been done to estimate the terminal height of hreballs (see sec¬ 
tion [^. In this paper, we have derived the terminal heights for MORP hreballs using 
newly developed as well as previously suggested methodology. This methodology had 


only been tested on several fully ablated hreballs with large /3 values by Gritsevich 
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and Popelenskaya (2008). We were particularly interested in determining whether 


this new mathematical approach works equally accurately with fully ablated hre- 
balls and meteorite-producing ones, and whether any additional modihcation could 
be applied to improve the accuracy and extend applicability range for the proposed 
methodology. 

The comparison of the derived terminal heights of MORP hreballs with their 
observed values brings us to the following conclusions. Using equation (21) it comes 
out that the accuracy shows a standard deviation of a = 1.52 km which seems to 
be high. However, a linear tendency is found. This last reason encouraged us to 
hnd a better adjustment. By applying newly suggested approximation in our study 
equation (24) we obtain terminal heights that £t very accurately the observed values; 
now the standard deviation for MORP data set turns out to be cr = 0.75 km. 

Besides, the excellent adjustment between derived and observed terminal height 
values justihes the methodology proposed here. We foresee a calculation of terminal 
height to be useful when the lower part of the trajectory was not instrumentally 
registered (e.g. limited held of view of the camera, shadowing of hreball trail by other 
objects, weather conditions, detection of meteoroids using radar techniques with 
limited observation beam etc.). It also brings critical knowledge into the problem 
when one needs to predict how long will be a total duration of the luminous hight 
or at which height a fireball produced by a meteoroid with given properties would 
terminate. Based on our investigations we can highly recommend the use of equation 
(24) also to solve inverse problem when terminal height and velocity are available 
from the observations, and parameters a and /3 need to be derived. 

Further observational and analytical research must be done to predict precisely 
the properties of any body along its atmospheric hight. 
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